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A method for detemiining a biomcchanical property of a feed, the method comprising the steps of: (a) subjecting the feed to infrared 
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Method for Determining Feed Quality 

This invention relates to a method for quantifying biomechanical properties of 
animal feed based on a correlation between the chemical and biomechanical 
properties of the feed, and to methods for objectively measuring the quality of 
5 animal feed, such as fodders including hay, pastures and forages. 

Diet is the major determinant of productivity of an animal. In the livestock 
industry, animals are farmed for meat, wool and other valuable products. The 
diet of farmed livestock is largely dictated by man and, given the effect of diet on 
animal production, it is highly desirable to optimise the diet of livestock to gain 
1 0 maximum benefit from the natural resource. 

Feed quality is one variable that has a major impact on animal productivity. In 

this respect, feed quality affects the amount of feed an animal will consume and 
_ the feeding value it gains from the feed consumed. In the case of cattle, sheep 

and other ruminants, feed quality depends on digestibTnty. chern^ 
1 5 (nutrient composition) and biomechanical attributes (namely how easy it is for an 

animal to chew the feed during ingestion and rumination). 

It is generally accepted that there are constraints on the intake of feed by 
ruminant animals, that the amount of useful energy obtained by a ruminant 
animal may fall short of the amount that the animal can potentially use, and that 

20 this would result in reduced productivity.. For example, the principal constraints 
to voluntary intake of fodders are resistance of fodder fibre to chewing and 
digestibility (provided that the intake is not otherwise constrained by low 
palatabllity, deleterious secondary compounds, or the inadequacy of essential 
nutrients). Differences between feeds, such as fodders, in their resistance to 

25 chewing are reflected in differences in biomechanical properties, including 
comminution energy, shear energy, compression energy, tensile strength, shear 
strength and intrinsic shear strength. 
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Hay is a common feed, and its quality is significantly affected by factors such as 
seasonal differences, haymaking practices and pasture composition. It has 
been shown In one recent survey that in some years as little as 11% of hay 
produced was good enough to promote liveweight gain in weaner sheep. This 
5 possibility of wide variation in measures of hay quality is a matter of increasing 
concern, and has given rise to a demand for a method of objective quality 
assessment. 

A hay quality system adopted in the United States of America uses a measure 
known as relative feed value (RFV) to distinguish between hays of different 
10 quality. The RFV is calculated from the dry matter digestibility, which is 
predicted from acid detergent fibre (ADF) content, and from the dry matter 
intake, which is predicted from neutral detergent fibre (NDF) content. 

The RFV based system suffers from a number of disadvantages. For examp/e. 

- - NDL^^^ by chemical methods 

15 which take several days to complete, and thus are expensive in terms of ' 
resources. 

While objective quality assessment and product specification has become an 
integral part of the production and marketing in domestic and export markets for 
the Australian grain, wool, meat and dairy industries, performance-based qualilty 
20 standards are not presently in place for feeds such as hays and other fodders. 
Consequently; 

(a) the feed buyer cannot be sure of getting value for money, and this bs 
likely to become increasingly important in respect of export markets if 
other exporting countries are able to guarantee standards for their 

25 product; 

(b) the feed producer cannot be sure of getting a higher price for a 
superior product; 
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(c) livestock producers are unable to objectively formulate rations or 
supplementary feeding regimes to achieve animal production targets; 
and 

(d) the market for animal feed tends to be unstable. 

5 Whilst the relationship between biomechanical properties of feed and feed 
quality is now accepted, there is a need for a convenient, inexpensive and 
relatively accurate assay method for feed to determine its quality. An accurate 
determination of feed quality allows for optimisation of feeding regimes and 
improved animal production for obvious economic gains. 

10 It is an object of this invention to overcome or at least partially alleviate the 
aforementioned problems and/or reduce the uncertainties and concomitant 
problems of the prior art systems for measuring the biomechanical properties of 
feed and hence determining feed quality. 



Thus, the present invention provides a method for determining a biomechanical 
1 5 property of a feed, the method comprising the steps of; 

(a) subjecting the feed to infrared radiation to obtain spectral data; 
and 

(b) using the spectral data to determine the biomechanical property; 

whereby the biomechanical property of the feed is determined on the basis of 
20 the bond energies of the chemical constituents of the feed. 

The spectral data may be used directly to determine the biomechanical property 
of the feed. Alternatively, the spectral data may be used to determine another, 
property of the feed and the other property is used to determine the 
biomechanical property on the basis of a correlation between the other property 
„25 andjhe biomechanical pro perty . 
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When the biomechanical property is determined via another property, the otiner 
property is preferably a chemical property of the feed such as the ADF content 
or the NDF content or the lignin content. 

There is a variety of biomechanical properties of the feed that may be 
5 determined. Preferably, the biomechanical properties are selected from the 
group comprising shear energy, compression energy, comminution energy, 
tensile strength, shear strength and intrinsic shear strength. 

The spectral data may comprise a reflectance spectrum at a combination of 
wavelengths or over a predetermined range of wavelengths such as 700nm- 
10 3000nm. or more preferably 1 100nm-2500nm. Preferably, the data obtained for 
the spectral range of 1850nm-1970nm is disregarded, this being the range over 
which water reflects strongly. 

The spectral data may be recorded at one or more wavelength intervals 
~ " throughout the spectral range. When the spectral data is a 
15 over a predetermined range it is preferably measured at 2nm intervals over the 

range. Of course, if so desired the spectral data may be measured at intervals 

other than 2nm. 

When the spectral data Is used to directly determine a biomechanical property, 
the biomechanical property is preferably determined by comparison of the 
20 spectral data with a calibration equation that reflects the relationship between 
reflectance and the biomechanical property. Preferably, the calibration curve Is 
determined on the basis of laboratory data establishing a correlation between 
reflectance and the biomechanical property. 

Thus, the present invention also provides a method for determining a 
25 biomechanical property of a feed, the method comprising the steps of; 

(a) subjecting the feed to infrared radiation to obtain spectral data; 

and — — '■ 
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(b) comparing the spectral data obtained in (a) with a calibration 
equation to determine the biomechanical property; 

whereby the biomechanical property of the feed is determined on the basis of 
the bond energies of the chemical constituents of the feed. 

5 The present invention also provides, a method for determining feed quality, the 
method comprising the steps of; 

(a) subjecting the feed tO: infrared radiation to obtain spectral data; 

(b) using the spectral data to determine a biomechanical property of 
the feed; and 

10 (c) using the value of the biomechanical property obtained in step (b) 

to determine feed quality; 

whereby the biomechanical property of the feed and thus the feed quality is 
determined on the basis of the bond energies of the chemical constituents of the 
feed. 

15 In one particular form, the method described immediately above may further 
comprise the determination of an additional property of the feed. The additional 
property may vary and preferably is selected from the group comprising the 
digestibility of the feed in vivo or in vitro, the ADF content or the NDF content, or 
the lignin content. 

20 The present invention is based on research establishing a strong correlation 
between the bond energies as they relate to the physical structure, and the 
biomechanical properties of feed. Once this correlation is established the bond 
energies of the chemical constituents, and in turn the biomechanical properties 
of the feed, can be determined using infrared spectroscopy. The biomechanical 
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properties quantified in this way are useful for accurately determining feed 
quality. 

In this respect, research resulting in the present invention has shown that the 
biomechanical attributes of feeds such as cereal and legume hays, straws, and 
mature, dry subterranean clovers are much more strongly related to animal 
performance than are digestibility or chemical composition of the feeds. 



Thus, comminution energy, the energy required to grind or comminute fodder 
material, has proved to be a very effective indicator of forage consumption 
constraint (FCC), which is the difference between the quantity of forage an 
10 animal should consume to satisfy its capacity to use energy (a theoretical 
maximum) and the actual voluntary dry matter intake achieved. 

Shear energy, the energy required to shear fodder material, and compression 
energy required to compress fodder material, are two biomechanical 



feed characters of fodders that are closely rerated to comminution energy and' 
15 which also are good predictors of FCC. 

In this respect, feed quality can be assessed in a number of ways. The forage 
consumption constraint (FCC) is one convenient measure of feed quality and 
equates to the difference between the quantity of the fodder that the animal 
would be attempting to consume to satisfy its capacity to use energy (theoretical 
20 maximum intake) and the voluntary forage consumption (VFC). 

Thus, the present invention also provides a method for determining feed quality, 
the method comprising the steps of; 

(a) subjecting the feed to infrared radiation to obtain spectral data; 

(b) using the spectral data to determine a biomechanical property of 
25 the feed; and 
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(c) using the value of the biomechanical property obtained in step (b) 
to deternrtine the forage consumption constraint (FCC) or voluntary 
feed consumption (VFC) as a measure of feed quality; 

whereby the biomechanical property of the feed and thus the feed quality is 
5 determined on the basis of the bond energies of the chemical constituents of the 
feed. 

The present invention is based on the finding that variations in biomechanical 
properties such as shear energy, comminution energy and compression energy 
are reflected in NIR spectra of fodders. This finding, together with recognition of 
10 the value of biomechanical characters for the prediction of FCC (and, in turn, the 
prediction of voluntary feed consumption (VFC) ) makes it possible for quicker, 
less expensive, more convenient and more reliable prediction of feed quality 
than hitherto known and predicted. 

Accbrdfngly. this'in^^ a ~ 

15 fodder, such as a forage, to meet a required animal performance; or (ii) 
predicting the VFC of a forage; or (iii) predicting the FCC of a forage, which 
method comprises subjecting a sample of the forage to NIR radiation and 
determining the reflectance at selected wavelengths. 

It has been found that the biomechanical properties, such as shear and 
2Q comminution energy values for a given fodder, correlate with the fodder's 
reflectance of infrared radiation. More specifically, the invention is based on 
research showing that: 

(a) NIR wavelengths at which reflectance (R), namely the second derivative of 
the logarithm of the inverse of R, correlates significantly with the variation in 
25 eriergy required to shear fodder materials are 1168nm. 1458nm, 1598nm, 

1718nm, 1828nm and 2048nm. For the prediction of fodder shear energy 
(y,. kJ.m'^) the following equation may be used: 
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= 19.95 + 10239.46 R.ies + 3623.49 R.^sb ' 4255.61 R.s^g - 
5319.88 Ri7is + 5148.38 Riege + 2452.05 R2048 

(b) NIR wavelengths at which the second derivative of the logarithm of the 
inverse of reflectance (R) correlates significantly with the variation in energy 

5 required to comminute fodder materials are 1138nm, 2018nm. 2128nm and 
2408nm. 

For the prediction of fodder comminute energy (yg. kJ.kg DM'^) the following 
equation is proposed: 

y2 = 231.42 + 18224.74 R^iag - 4955.12 Rgoia - 3005.37 R2128 + 
10 4290.18 R2408 

(c) NIR wavelengths at which the second derivative of the logarithm of the 
inverse of reflectance (R) correlates significantly with the variation in 
compression energy are 1 268hm, n 588nm7~ 1 728nm - 
prediction of compression energy (yg. kJ.kgDM'^) the following equation may 

15 be used: 

y3 = -0,71 - 91 1 ,04 R1268 + 1 12,57 R^^^ - 79.48 R1728 - 28.02 R227B 

(d) NIR wavelengths at which the second derivative of the logarithm of the 
inverse of reflectance (R) correlates significantly with variation in in vivo 
digestibility of dry matter (DMD) (y4, %) is 1158nm, 1238nm, 1668nm, 

20 1908nm, 1918nm. and 2248nm. For prediction of the DMD (y4. %) of a 

fodder the following equation is proposed: 

y4 = 46.62 + 8162.72 R^ge - 8799,69 R1238 + 1249.01 Riges + 519.46 
R1908 - 367,08 R1918 - 161 ,84 R^^^ 

(e) NIR wavelength at which the second derivative of the logarithm of the 
-25 inverse of reflectance--(R)--correlates„significantly_withLjyar^^^ 
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digestibility of dry matter (IVDMD) is 1698nm. 1748nm. 1908nm, 1918nm 
and 2158nm. For prediction of the DMD {in vitro) of a fodder the following 
equation is proposed: 

Ys = 63.43 - 2186.89 8^698 ' 1491.99 R1748 + 981.30 R^gos - 556.01 
5 R1918 + 2003.05 R2158 

Accordingly, in a preferred method according to this invention, the infrared 
wavelengths at which reflectance is measured comprise one or more of the 
following: 1168nm, 1458nm, 1598nm, 1718nm. 1828nm, 2048nm, 1138nm. 
2018nm, 2128nm. 2408nm, 1268nm. 1588nm, 1728nm. 2278nm, 1158nm. 
10 1238nm. 1668nm, 1908nm. 224Bnm. 1698nm. 1748nm, 1918nm and 

2158nm. 

It will be understood that the foregoing are wavelengths at which the 
strongest correlations have been observed, and the possibility of useful 
correlations being obsefve'd a^^ 

15 Essentially, it can be shown that in the same way that a decrease in 
comminution energy is reflected by a decrease in forage consumption 
constraint, there is also a linear relationship between comminution energy or 
shear energy and the consumption constraint of a fodder. Thus, the use of 
NIR spectra, in conjunction with the equations detailed at paragraphs (a) to 

20 (e) above, permits estimation of the VFC of a fodder, which together with 
estimates of digestibility (conveniently obtained from NIR spectra) can be 
expected to provide a valuable basis for performance-based quality 
standards for fodders. 

It is to be appreciated that the intention of this invention is to offer a quick, 
25 reliable and relatively inexpensive means of obtaining information -from 
which the fodder producer and user, such as purchaser, might make 
informed judgements about the market value of a given fodder sample 
l-'ilat'ive'to altem^^^^ of its suitability for a-particular-purpose. 
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Conceivably. fodder quality predictions obtained by the method of this 
invention could be a useful component of. or used in conjunction with, for 
example. Decision Support Software (DSS) packages designed to assist 
livestock management. 

i It is further envisaged that by combining NIR measurements made by a 
remote sensing system, such as Undsat, with data from a Geographical 
Information System, the invention will provide a means of making reliable 
predictions of pasture quality. These predictions, together with predictions 
of feed intake and animal performance, should then provide a useful basis 
for strategies of supplementary feeding to improve performance in grazing 
ruminants. 



The present invention also provides for a spectrometer configured to 
determine biomechanical properties and/or quality of feed according to the 

- methods _oltJie_presennnyention.^^^ ^^^^^ the spectrorrieter includ_es a 

15 data processing means which enables the spectrometer to receive a feed 
sample and quantify either or both the biomechanical properties of the feed 
and the quality of the feed. In one particular form the data processing 
means includes a calibration equation to facilitate the determination of the 
feed quality or biomechanical property. 

20 The invention will now be described with reference to the following 
examples. The description of the examples is in no way to limit the 
generality of the preceding paragraphs. 
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EXAMPLES 

The energy of molecular vibrations correspond to the energy of the infrared 
spectrum of the electromagnetic spectrum, and these molecular vibrations may 
be detected and measured in the wavelength range of the infrared spectrum. 
5 Functional groups in molecules have vibration frequencies that are characteristic 
of that functional group and that are within well-defined regions of the infrared 
spectrum. 

For organic compounds the principal analytical features of the near infrared 
10 (NIR) spectrum are due to absorbance of radiant energy by bonds between 
hydrogen, carbon, nitrogen, oxygen or with sulphur, phosphorus and metal 
halides. When organic compounds are irradiated with infrared radiation at 
wavelengths between 700 and 3000nm part of the' incident radiation is absorbed 
and the remainder is reflected, refracted or transmitted by the sample. Most 
15 -quantitative reflectance. analyses. are. mM^.^^i^^^^^ w^^velength range of 11 00 to 
2600nm. The amount of energy absorbed or diffusely reflected at any given 
wavelength in this wavelength range is related to the chemical composition of 
the organic compound. NIR spectroscopy uses detectors to measure the 
amount of radiation that is diffusely reflected by the irradiated sample. 

20 

NIR spectroscopic analysis is an analytical procedure calibrated to a primary 
reference method. Calibration in NIR spectroscopy (NIRS) relies on similarities 
among the spectra, and analytical properties of interest in the reference 
samples. In this example the analytical properties of interest were the 
25 biomechanical characters of forages, and the procedure that was adopted in this 
example was as follows: 

a) prediction of biomechanical characters of a range of grasses using NIR 
spectroscopy was established by developing a calibration equation(s) from 
30 laboratory determined values of a set of reference samples. 
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b) validation of the equation(s) either by using laboratory determined values of a 
separate set of samples, or by a cross-validation procedure using the 
laboratory determined values of the reference samples. 

c) using the NIRS-predicted values for biomechanical characters of the forages 
5 and for digestibility of the forages, forage consumption constraint (FCC) was 

predicted, and in turn voluntary feed consumption (VFC) was predicted. 

d) the predicted FCC and VFC were compared with actual data from groups of 
animals fed each of these forages. 

10 Example A: Developing a calibration equation to predict biomechanical 
properties of herbage: 

The samples used in this example were a range of varieties of Panicum spp. 
harvested at a range of plant maturities throughout the growing season (Tabfle 
15 1). Each of the samples was dried and chaffed, and then fed to groups of sheep 
(8 sheep per group) which were penned individually, to determine 7/7- wVo dry 
matter digestibility (DMD), VFC and FCC. Samples of the hays were stored for 
laboratory analyses. 

20 Biomechanical properties of the forages were determined using published 
methods: the energies required to shear or compress the forages according to 
Baker, Klein, de Boer and Purser (Genotypes of dry, mature subterranean clover 
differ in shear energy. Proceedings of the XVII International Grassland 
Congress 1 993. pp 592-593.) and the energy required to comminute the forages 

25 according to Weston and Davis (The significance of four forage characters as 
constraints to voluntary intake. Proceedings of the Third International 
Symposium on the Nutrition of Herbivores, Penang, 1991). In vitro digestibility of 
dry matter (IVDMD) was determined by the pepsin-cellulase technique as 
modified by Klein and Baker (Composition of the fractions of dry, mature 

30 subterranean clover digested in vivo and in vitro. Proceedings of the XVII 
International Grassland Congress 1993. pp593-595.). 
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There are several ways to process samples for NIRS analysis, and in this 
example the samples were ground through a cyclone mill with a 1 mm screen 
and equilibrated at 25°C for at least 24h before NIRS analysis. The samples 
were scanned by a monochromating near infrared reflectance 
5 spectrophotometer (Perstorp NIRS 6500) and the absorption spectra recorded 
for the range 1100 to 2500nm at 2nm intervals. The spectral range 1850 to 
1970nm, where water absorbs strongly, was disregarded in further analysis of 
the spectra! data. 

10 For NIRS analysis the samples were divided into two groups: one group to be 
used as a 'calibration* set to establish a prediction equation, and a second 
group, the Validation' set. to be used to validate the prediction equation. There 
are a number of ways to select the samples for each set. In this example the 
samples were ranked according to each of the characters that were to be 

15 predicted and every other sample was selected for the calibration set (33 
samples) and the validation set (32 samples).— Thus, for each characterJhat. was _ 
evaluated, a different selection was made from the 65 samples to establish the 
respective calibration and validation sample sets. 

20 The ranges, mean, median and variation in the laboratory-determined values for 
each of the characters of interest in the calibration and validation sets are listed 
in Table 2. 

The software for scanning, mathematical processing and statistical analysis were 
25 supplied with the spectrophotometer by the manufacturers. The spectral data 
were transformed by taking the second derivative of the logarithm of the inverse 
of the reflectance (R) at each wavelength (d" log (1/R)). The similarities amongst 
the spectra (Figure 1) of the samples in the validation and calibration sets were 
determined using principal components scores to rank the spectra according to 
30 the Mahalanobis distance from the average of the spectra. The Mahalanobis 
distance values were standardized by dividing them by their average value, and 
gveTellenote 
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Calibration equations were developed using the calibration samples by 
regressing the data from the laboratory analyses of each biomechanical property 
against the corresponding transformed spectral data using the following 
5 mathematical methods: 

a) Stepwise linear regression 

b) Step-up linear regression 

c) Principal components regression (PGR) 

10 d) Partial least squares regression (PLS), and 

e) Modified partial least squares regression (MPLS). 

Stepwise calibrations were developed for each calibration set of samples using 
the mathematical treatments of the spectral data 2.2.2; 2,5,5; 2.10,5; and 

15 2,10.10; where the first number denotes that the second derivative was used, 

the second Indicates that second derivatives of the spectral data (determined at 

2nm Intervals) were taken at intervals of 4, 10 or 20nm, and the third indicates 
that the function was smoothed using the 'boxcar* method over Intervals of 
wavelength of 4. 10. or 20nm (Table 4a). Likewise step-up calibrations were 

20 developed for each calibration set with up to 6 terms in each calibration equation 
using mathematical treatments 2,2,2; 2.5.5; 2,10,5; and 2,10.10 CTable 4b). 
Calibrations developed for each calibration set using principal components 
regression, partial least squares regression, or modified partial least squares 
regression each were developed using mathematical treatments 2,5.5 and 

25 2.10.10 CTable 4c). 

In developing the calibration equations In the stepwise and step-up regressions, 
only wavelengths with partial F-statistIc of more than 8 were accepted for the 
models. 



30 



For each calibration using each calibration set the following calibration statistics 
"wisre'deterrnined: 
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a) Squared multiple correlation coefficient (R^), an indication of the proportion of 
the variation in the calibration set that is adequately modelled by the 
calibration equation. 

b) The standard error of calibration (SEC) together with its confidence interval (± 
5 CL). which is the standard deviation for the residuals due to difference 

between the laboratory determined (reference) and the NIR predicted values 
for samples within the calibration set 

Once the calibration equations were developed, each equation was validated by 
10 using it to predict the respective biomechanlcal property values for each sample 
in the validation sample set. For each calibration equation the following 
validation statistics were determined: 

a) Simple linear correlation coefficient (r^) between the laboratory determined 
and NIR predicted values. 
15 b) The bias (or systematic error) in the regression relationship between the 
— - laboratory determined (reference) and fsJIR predicted values. 

c) The confidence limits of the bias in the regression relationship between the 
laboratory determined (reference) and NIR predicted values. 

d) The standard error of prediction, corrected for bias (SEP(C)). which 
20 represents the unexplained error of the prediction, the deviation of the 

differences between laboratory determined and NIR predicted values, 

e) The coefficient of determination, or slope (p), and y-intercept (a) of the linear 
regression relationship between the laboratory determined and NIR predicted 
values. 

25 f) The residua! standard deviation (RSD) of the linear regression relationship 
between the laboratory determined and NIR predicted values. 

In addition, the calibration equations were validated using a procedures of cross- 
validation. These are procedures where every sample in the calibration set was 

30 used once for prediction, and the standard error of validation corrected for bias 
(SE\/(C), for stepwise and step-up regressions) and cross-validation (SECV, for 

multivariate regressions)-Gan be determined 
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Calibration equations for each biomechanical character were selected using the 
following criteria: 

a) Lowest partial F-ratio. highest lowest SEC and. for PCR. PLS and MPLS, 
lowest SEV(C) (or. for multivariate regressions, SECV) 

b) Highest r^. lowest bias and I bias I < bias confidence limit, lowest SEP(C), p 
closest to 1.0, a closest to 0. and lowest RSD. As well. SEP(C) was 
compared with the standard error of laboratory determined values amongst 
all 65 samples, listed in Table 5. 

Calibration equations were similarly established to predict in vivo digestibility and 
in vitro digestibility. The coefficients for each wavelength in the selected 
calibration equations from stepwise or step-up regression analyses are listed in 
Table 6a, and those from multivariate analyses are listed in Table 6b. 



S'"^P'® linear correlation coefficient (r^) between the laboratory determined and 
NIR predicted values for each of the biomechanical characters (energies 
required to shear, comminute or compress) and digestibility of dry matter 
determined in vivo or in vitro of the samples in the validation set are shown in 
20 Figures 2a. 2b. and 2c. The NIR predicted values are predicted using calibration 
equations that best met the criteria listed above. 

Example B: Prediction of FCC and VFC using NIR determinations of energy 
required to shear and in vivo digestibility: 

25 

To demonstrate the prediction of voluntary feed consumption using NIR 
determined values for a biomechanical character and digestibility of forages, 
samples of Panicum spp. hay were selected which were common to both of the 
validation sample sets used to establish the NIR prediction equations for energy 
30 required to shear and in vivo digestibility. The hays represented the range o-f 
varieties in the sample set, and are listed in Table 7. The samples were 
scanned ~ by" the— same-spectrophotometer - that-was~crsed-to— establish-the" 
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calibration equations, and the absorption spectra were recorded in the range 
1 1 00 to 250Gnm at 2nm intervals. Values for energy required to shear and in 
vivo digestibility were predicted from calibration equations (Tables 4a, 4b and 
4c) using the recorded spectral data. 

5 

These values then were used to estimate FCC from the relationship between 
biomechanical character(s) and FCC of the range of forages used by Weston 
and Davis (1991). Energy required to shear the forages used by Weston and 
Davis was determined according to Baker et a/. (1993). The relationship 
1 0 between the energy required to shear these forages (kJ/m^) and FCC (g organic 
matter (OM) / d / kg metabolic body weight (MBW)) was described by the 
relationship: 

Energy required to shear (x) = -26.1 3 + 5.53 (FCC (y)) 
where R = 0.92; RSD = 8.70; N = 1 3; P < 0.0001 . 

FCC from this relationship, and /n j^/vo digestibility^^^ 

used to estimate VFC. as the difference between the animal's capacity to use 
energy (as defined by Weston and Davis, 1991) and FCC. These data are 
summarised in Table 8. 

20 

VFC predicted in this way explained most of the variation in actual VFC (R = 
0.87; RSD = 5.04; P = 0.023) (Figure 3). 
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Table 2. Summary statistics for each calibration and validation set 





Energy required 
to shear 

(kJ/m^ 


Energy required 
to comminute 
fkJ/kq DM) 


Energy required 
to compress 
(kJ/kq DM) 


OigestibitJty of 
dry matter ^n vivo 
f%) 


Digestibility of dry " 
matter in vitro 


ICalibratton sot 





mean 


15.48 


134.9 


1 median 


15.17 


133.8 


I maximum 


20.95 


216.5 


minimum 


10.80 


72.5 


j standard deviation 


2.572 


37.50 


Valldatfon set 






mean 


15.43 


130.9 


median 


15.20 


126.3 


maxtmum 


20.43 


205.2 


minimum 


10.94 


54.5 


standard deviation 


2.444 


37.50 



ICalibratlon set 
mean 
median 
maximum 
minimum 

standard deviation 
iValtdatlon set 
mean 
median 
maximum 
minimum 

standard deviation 



Energy fequrre<t to comminute 



3,70 


55.7 


53.3 


3.65 


56.0 


55.1 


4,39 


64.0 


63.0 


3.25 


43.0 


39.8 


0.265 


5.73 


6.37 


3.78 * 


55.6 


52.7 


3.75 


56.5 


53.3 


4.24 


64.0 


63.0 


3.34 


47.0 


40.1 


0.229 


5.36 


7.01 



15.01 
14.76 
19.97 
10.80 
2.444 


133.1 
129.5 
216.5 
54.5 
38.62 


3.69 
3.70 
4.18 
3.25 
0.227 


55.7 
57.0 
64.0 
43.0 
5.64 


52.8 
54.7 
63.0 
398 
7.06 


15.92 
15.97 
20.95 
11.46 
2.490 


132.9 
130.2 
205.2 
60.7 
36.20 


3.79 
3,79 
4.39 
3.34 
0.263 


55.6 
55.5 
64.0 
47.0 
5.47 


53.2 
54.7 
62.5 
40.9 
6.92 



jCallbratlon set 

mean 

median 

maximum 

minimum 

standard deviation 
{Validation set 
mean 
median 
maximum 
minimum 

standard deviation 



Enerpy repulred lo compress 



15.28 
15.07 
19.97 
10.80 
2.477 

15.64 
15.42 
20,95 
11.46 
2.530 



128.1- 
126.4 
204,0 
54.5 
38.00 

138.0 
132.5 
216.5 
60.7 
36.39 



3.74 
3.72 
4.39 
3.25 
0.261 

3.74 
3.72 
4.24 
3.34 
0.240 



-56.3 
57.0 
64,0 
47.0 
5.15 

55.0 
54.5 
64.0 
43.0 
5.87 



53.8 
54.7 
63.0 
39.8 
6.64 

52.2 
54.5 
62.0 
40.1 
7.26 



Dlgeet<bM>ty of dry matter in vivo 



(Calibration set 
mean 
median 
maximum 
minimum 

standard deviation 
jValidatlon set 
mean 
median 
maximum 
minimum 

standard deviation 



15.14 


133.9 


3.74 


55.5 


53.1 


15.17 


128.4 


3.72 


56.0 


55.1 


20.95 


216.5 


4.24 


64.0 


63.0 


10.80 


60.7 


3.25 


43.0 


40.1 


2.528 


36.39 


0.247 


5.73 


7.33 


15.78 


132.1 


3.75 


55.7 


52.9 


15.20 


134.7 


3.72 


56.5 


54.4 


20.37 


205.2 


4.39 


64.0 


63.0 


10.94 


54.5 


3.34 


47.0 


39.8 


2.446 


38.70 


0.255 


5.36 


6.63 
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Table 2 (cont'd), 
set 



Summary statistics for each calibration and validation 





Energy required to 


Energy required 


Energy required 


Digestibility of 


Digestibility of dry 




shear 


to comminute 


to compress 


dry matter in 


matter in vitro 




(IcJ/m^) 


(kJ/kg DM) 


(kj/kg DM) 


vivo 


(%) 








(%) 




^.t:.^ ■..•.DIoefit1b«ltyt)t:drv^iTiim«rr^^^ 


Calibration aat 












mean 


14.70 


131.3 


3.75 


55.6 


53.0 


median 


14.34 


129.3 


3.71 


56.0 


54.7 


maximum 


19.36 


216.5 


4.39 


64.0 


63.0 


minimum 


10.80 


54.5 


3.25 


43.0 


40.1 


standard 


2.235 


42.58 


0.241 


5.78 


6.94 


deviation 












Validation sat 










53.0 


mean 


16.19 


134.6 


3.74 


55.6 


median 


16.16 


133.8 


3.74 


56.0 


54.7 


maximum 


20.95 


194.8 


4.24 


64.0 


63.0 


minimum 


10.94 


65.7 


3.34 


47.0 


39.8 


standard 


2.538 


31.64 


0.260 


5.33 


7.05 


deviation 













Table 3. 



Mahalanobis distances 



Mean 



— Median 



Range 



For full sample set: 

For calibration sets for: 
Energy required to shear 
Energy required to comminute 
Energy required to compress 
Digestibility of dry matter in vivo 
Digestibility of dry matter in vitro 



0.655 

0.588 
0.718 
0.757 
0.673 
0.645 



0.623 

0.549 
0.676 
0.760 
0.634 
0.574 



0.203 - 1.983 

0.171 - 1.646 
0.350- 1.553 
0.188- 1.440 
0.389 - 1 .547 
0.185- 1.178 
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Table 4a. Calibration and validation statistics 



PCT/AU96/00776 



energy re<ru>red to shear 



2.2,2 



Stepwis, Regression 



2.5.5 



Lowest partial F-ratio 

SEC 
SEC CL 
SEV(C) 

Bias 
Bias CL 
SEP (C) 
Slope 
Intercept 
R.S.D. 

IBiasI • Bias CL 
jBiasI < Bias CL? 



2.10.5 2.10.10 



10.27 
0.798 
1.155 
1.493 
1.230 
0.368 
0.690 
1.484 
1.500 
0.604 
6.340 
1.627 
-0.794 
Yes 



6.18 
0.787 
1.188 
1.535 
1,306 
0.625 
0.710 
1.527 
1.540 
0.617 
5.440 
1.627 
-0.817 
Yes 



8.27 
0.795 
1.166 
1.507 
1.273 
0.520 
0.700 
1.498 
1.520 
0.598 
5.640 
1.484 
-0.798 
Yes 



4.70 
0.780 
1.207 
1.560 
1.322 
0.495 
0.720 
1.551 
1.570 
0.758 
3.710 
1.476 
-0.831 
Yes 



\number of terms 



Stepwise Regression 





2.2.2 


2.5,5 


2.10,5 


2.10.10 


Lowest partial F-ratro 


5.54 


4.45 


16.55 


10.89 


R^ 


0.910 


0.602 


0.818 


0.831 


SEC 


1 1 .626 


17.281 


16.546 


15.980 


SEC CL 


1.493 


1.535 


1.507 


1.560 


SEV(C) 


13.103 


18.040 


.. 17.587 . 


17.100 




0.363 


0.429 


0.374 


0.213 


Bias 


6.980 


10.370 


9.930 


9.590 


Bias CL 


14.941 


22.209 


21.264 


20.537 


SEP (C) 


15.110 


22.460 


21.510 


20.770 


Slope 


0.530 


0.575 


0.607 


0.417 


Intercept 


58.300 


48.900 


48.600 


74.600 


R.S.D. 


28.900 


27.360 


28.650 


32.120 


IBiasI - Bias CL 


-7.961 


-11.839 


-11.334 


-10.947 


IBiasj < Bias CL? 


Yes 


Yes 


Yes 


Yes 



{number of terms 



Stepwise Regression 





2,2.2 


2.5.5 


2.10.5 


2.10.10 


1 Lowest partial F-ratio 


5.05 


4.44 


7.90 


16.19 


R^ 


0.784 


0.500 


0.525 


0.534 


SEC 


0.121 


0.209 


0.204 


0.202 


SECCL 


1.493 


1.535 


1,507 


1.560 


SEV{C) 


0.135 


0.224 


0.217 


0.215 




0.069 


0.113 


0.008 


0.067 


1 Bias 


0.070 


0.130 


0.120 


0.120 


Bias CL 


0.156 


0.269 


0.262 


0.260 


SEP (C) 


0.160 


0.270 


0.270 


0,260 


Slope 


0.180 


-0.080 


0.314 


0.21 1 


Intercept 


3.060 


4.030 


2.580 


2.960 


R.S.D. 


0.229 


0.229 


0.227 


0.232 


IBiasI - Bias CL 


-0.086 


-0.139 


-0.142 


-0.140 


IBiasI < Bias CL? 


Yes 


Yes 


Yes 


Yes 


{number of terms 


6 


4 


4 


4 1 
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Table 4a (cont'd) 



; DIa stiblltty otdrjr matter ^Wvo ^-^ - ;^ 






Stepwis 


Regression 






2.2.2 


2.5.5 


2.10.5 


2.10.10 


Lowest partial F-ratio 


7.63 


20.68 


4.28 


6.08 




0.934 


0.917 


0.914 


0.921 


SEC 


1.107 


1.236 


1.258 


1.207 


SEC CL 


1.493 


1.535 


1.507 


1.560 


SEV{C) 


1.215 


1.368 


1.341 


1.284 




0.654 


0.881 


0.878 


0.876 


Bias 


1.070 


0.890 


0,910 


0.900 


Bias CL 


0.156 


0.269 


0.262 


0.260 


SEP (C) 


2.320 


1.940 


1.980 


1.960 


Slope 


0.705 


0.878 


0.840 


0.827 


Intercept 


16.500 


6.690 


6.640 


9.340 


R.S.D. 


3.153 


1.852 


1.873 


1.888 


(Biasl - Bias CL 


0.914 


0,621 


0.648 


0.640 


IBiasI < Bias CL? 


No 


No 


No 


No 


number of terms 


6 


6 


6 


5 










Stepwise Regression 






2.2.2.1 


2,5.5 


2.10.5 


2,10.10 


Lowest partial F-ratio 


7.68 


11.84 


4.33 


6.31 


R= 


0.935 


0.933 


0.915 


0.922 


SEC 


1.808 


1.751 


2.052 


1,974 


SEC CL 


1.493 


1.535 


1.507 


1.560 


SEy(C) 


1.984 


1.981 


2.186 


2.100 




0.699 " 


" 0.847 - 


0.743 


0.736 


Bias 


1.080 


1.050 


1.230 


1.180 


Bias CL 


2,324 


2.250 


2.637 


2.537 


SEP (C) 


2,340 


2.280 


2.670 


2.570 


Slope 


0.839 


0.962 


0.775 


0.763 


Intercept 


8.790 


1.650 


12.200 


12.700 


R.S.D. 


3.805 


3.794 


3.805 


2.719 


|Bias| - Bias CL 


-1.244 


-1.200 


-1.407 


-1,357 


IBiasI < Bias CL? 


Yes 


Yes 


Yes 


Yes 


number of terms 


6 


5 


6 


5 
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Table 4c. Calibration and validation statistics (multivariate 



regressions) 



• ^>:Eher.9V-reqUlred' tb'shear- "^^^^^ 




pnR PLS MPLS 


2.5,5 2.10.10 


2.5.5 2.10.10 


2.5.5 2,10.10 


SECV 

Bias 
Bias CL 
SEP (C) 
Slope 
Intercept 
R.S.D. 

IBiasI - Bias CL 
jBiasI < Bias CL? 


0.847 0.752 
1 .036 1 .290 
1.199 1.493 
1.750 1.592 

0.5441 0.4876 
0.620 0.770 
1.331 1.658 
1.350 1.680 

0.6540 0.6850 

5.2900 4.8100 
1.671 1.776 

.0.71 1 -0.888 
Yes Yes 


0.639 0.601 
1 .545 1 .624 
1.788 1-879 
1 .788 1 .933 

0.4938 0.4157 
0.930 0.970 
1 .986 2.087 
2.010 2.110 

0.7390 0.6270 

3.7500 5.4500 
1 .761 1 .892 

-1.056 -1.117 
Yes Yes 


0.601 0.582 
1,550 1.586 
1.793 1.835 
1 .600 1 .583 

0.3080 0.3563 
0,930 0.950 
1.992 2.038 
2.020 2.060 

0.5220 0.6000 

7.3100 6.0200 
2.065 1 .992 

-1.062 -1.088 
Yes Yes 








PGR PLS MPLS 


2,5.5 2.10,10 


2.5,5 2.10.10 


2.5.5 2.10.10 


R^ 
SEC 
SEC CL 
SECV 

Bias 

-BiasCL. 
SEP (C) 
Slope 
Intercept 
R.S.D. 

IBiasI • Bias CL 
IBiasI < Bras CL? 


0.584 0.574 
23.378 23.682 
27.048 27.400 
26-030 26.121 
0.349 0.337 
14.030 14.210 
.__30.044 30.435 
"30!390 " 30.790 " 
0.649 0.636 
37.3 39.2 
28.246 28.676 
-16.014 -16.225 
Yes Yes 


0.605 0.595 
22.788 23.075 
26.366 26.698 
25.548 25.683 

0.332 0.325 
13.670 13.840 
29.286 29.655 
" 29-620 30.000" 

0.644 0.632 

38.1 39.7 
28.714 28.884 
-15.616 -15-815 
Yes Yes 


0.556 0.558 
24.164 24.101 
27,958 27.885 

26.409 26.252 
0.33 0.328 

14.500 14.460 
31-055 30.974 

31.410 -31.330 
0.657 0.638 
36.6 39.1 

28.651 28.900 
-16.555 -16.514 
Yes Yes 




-;V:;;^^;:Eneft}Y reiqulred4o:compress*^ • ^.^ ^ ^ --r^^: ^.ym^^y^^ 




PGR PLS MPLS 


2.5.5 2.10,10 


2.5.5 2.10.10 


2.5.5 2,10.10 


R^ 
SEC 
SEC CL 
SECV 

Bias 
Bias CL 
SEP (C) 
Slope 
Intercept 
R.S.D. 

IBiasI ' Bias CL 
IBiasI < Bias CL? 


0.251 0.160 
0.225 0.241 
0.260 0.279 
0.299 0.277 
0.0220 0.0120 
0.140 0.140 
0.289 0.310 
0.290 0.310 
0.2290 0.2270 
2.8900 2.9000 
0.235 0.236 
-0.149 -0.170 
Yes Yes 


0.231 0.208 
0.265 0.269 
0.307 0.311 
0.301 0,307 
0.0130 . 0.0090 
0.160 0.160 
0.341 0.346 
0.340 . 0.350 
0.1530 0.1330 
3.1700 3.2500 
0.236 0.237 
-0.181 -0.186 
1 Yes Yes 


0.038 0.040 
0.260 0.260 
0.301 0.301 
0.301 0.299 
0.0060 0.0080 
0.160 0.160 
0.334 0.334 
0,340 0.340 
0.2290 0.2590 
2.8900 2.7700 
0.237 0.237 
-0.174 -0.174 
Yes Yes 
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Table 4c (cont'd) Calibration and validation statistics 
(multi variate regressions) 



_PCR pLS 
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Table 5. 



-Si- 
Standard error of laboratory determination (SEL) 





Energy 


Energy 


Energy 


Olgestlbitity 




required to 


required to 


required to 


of dry 




shear 


comminute 


compress 


matter 










in vivo 




fkJ/m*) 


(kJ/kg DM) 


(kJ/kg OM) 


f%) 


Mean SEL (n^5) 


0.796 


5.630 


0.078 


not svmdkUa 


Median SEL 


0.788 


5.492 


0.085 


notevstlsM 


Maximum SQL 


Z044 


13.098 


0.211 




Minimum SEL 


0.114 


0.750 


0.019 


notavsikbie 


SEL CL (using mean SEL) 


1.035 


7.319 


0.101 


notmdable 


SEL CL (using median SEL) 


1.024 


7.140 


0.111 





DIgestlbllrty 
of dry 
matter 

in vitro 

J2y_ 



0.314 
0.270 
1.126 
0.005 
0.408 
0.351 
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Tabie 6a. Components of possible prediction equations 
froRi stepwise and step-up rsgrossion analyses. 







Ccofilcfant Wavofcmctti 






ftrcd to shear 


Mattuimatteai trocCrrtcjrt 


2 7,7 fg tcjTTO) 


Step-up 
Z5.5 a terms) 




-'S^udl 1S58 
3323.0 14a 
-S3ig^ 1718 

9i4a38 ie2s 

102^.43 11^ 


28.09 

1033.77 2048 
700.12 1638 


Eiujygy rcquirod to eompross 




ZtO.IOMtcjmo) 


Stop-up 
ZlO,5f3tcnTO) 




^71 

•38102 2278 
11^57 13W 
-79.48 1T28 
-911.04 12OT 


Z48 

-31 -OS 1728 
-108.89 1546 
-4<a^ 12^ 


Enomv require to comminutQ 


Rogrosston qnstpsLa 
MoitKamoClcal troatmonft 


Z10.5 (4tcrTTO) 


Stop-up 
Z10.5 (1 tOTi) 




23lT42 

•3008137 2128 
4Sai9 24{» 
.4SS3.12 2»8 
18224.74 1138 


^.08 

-1321;33 2238 


DSgostlbliny of dry msrtter /rr v#t*no 




25.5 (6 XisrmQ) 


Stop-up 
Z10.S Ptsms) 




4a62 

-337.08 1918 
-87S9.68 1238 
8182-72 1158 

12^.01 ^em 

319.48 1OT 
-181.04 2248 


48.18 

-61Z43 1696 
SZa2 1418 
.S43.77 1618 


tMqostfblltty of dry matter 2r? vftra 


Regrossksn analysis 
Mattvamctlcal tro^rtmsnt 


25.5 (5 


Step-up 
Z10.10(4temi3) 




83.43 

-558.01 1918 

CT1JS0 isoe 

-21M^ 1688 
2003^06 2158 
-1491^ 1748 


54.29 

-1171.70 1698 
311.12 1418 
-2857.69 1818 
-2319.81 1228 
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THE CLAIMS of the invention are as follows: 

1. A method for determining a biomechanical property of a feed, the method 
comprising the steps of: 

(a) subjecting the feed to infrared radiation to obtain spectral data; 
5 and 

(b) using the spectral data to determine the biomechanical property; 

whereby, the biomechanical property of the feed is determined on the basis 
of the bond energies of the chemical constituents of the feed. 

2. A method according to claim 1 wherein the biomechanical property of the 
1 0 feed is determined directly from the spectral data. 

3. A" 'method according to clatrri 1 "wherein" the 'speictral data is used to 
determine another property of the feed and the other property is used to 
determine the biomechanical property on the basis of a correlation between 
the other property and the biomechanical property. 

15 4. A method according to claim 3 wherein the other property is ADF content, 
NDF content or lignin content. 

5. A method according to claim 1 or claim 2 wherein the spectral data is a 
reflectance spectrum over a predetermined range of wavelengths. 

6. A method according to claim 5 wherein the predetermined range is 
20 approximately 700nm to 3000nm, 

7. A method according to claim 5 wherein the predetermined range is 
approximately 1 1 0Onm to 2500nm. 
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8. A method according to any one of claims 5 to 7 wherein the data obtained 
for the spectral range of approximately 1 850nm to 1 970nm is disregarded. 

9. A method according to any one of claims 5 to 8 wherein the spectral data is 
recorded at 2nm intervals over the predetermined range. 

5 1 0. A method according to claim 1 or claim 2 wherein the reflectance reading is 
taken at a combination of wavelengths. 

11. A method according to claim 1 0 wherein the combination of wavelengths is 
selected from the group comprising: 1168nm. 1458nm, 1598nm. 1718nm, 
1828nm, 2048nm. 1138nm, 2018nm, 2128nm. 2408nm, 1268nm. 1588nm, 

10 1728nm. 2278nm. 1158nm. 1238nm. 1668nm, 1908nm, 2248nm, 1698nm, 

1748nm, 1918nm and 2158nm. 

12. A method according to claim 10 wherein the combination of wavelengths is 
11 68nm, , 1 458nm, „.1 598nm,. 1 71 8nm, . 1 828nm. ..and . 2048nm ._and the 
biomechanical property is shear energy. 

15 13. A method according to claim 10 wherein the combination of wavelengths is 
1268nm. 1588nm. 1728nm and 2278nm and the biomechanical property is 
compression energy. 

A method according to claim 1 0 wherein the combination of wavelengths is 
1138nm. 201 8nm. 2128nm and 2408nm and the biomechanical property is 
comminution energy. 

15. A method for determining a biomechanical property of a feed, the method 
comprising the steps of: 

subjecting the feed to infrared radiation to obtain spectral data; 
and 



14. 

20 



(a) 
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(b) comparing the spectral data obtained in (a) with a calibration 
equation to determine the biomechanical property; 

whereby, the biomechanical property of the feed is determined on the basis 
of the bond energies of the chemical constituents of the feed. 

5 16, A method according to claim 15 wherein the calibration equation is y, = 
19.95 + 10239.46 + 3623.49 R^^^b - 4255.61 Rigge - 5319.88 R^y^^ + 

5148.38 Ri82B + 2452.05 R2046 the biomechanical property is shear 
energy(y,). 

17. A method according to claim 15 wherein the calibration equation is yg = 
10 231.42 + 18224.74 R^^^s • 4955.12 Rgoie " 3005.37 R212B + 4290.18 R24oa 

and the biomechanical property is comminution energy (yg). 

18. A method according to claim 15 wherein the calibration equation is ys = - 
0.71- 911.04 Ri26e +112.57 R,5e8 - 79,48 Ri72B -- 28.02 _R2278 and the 

biomechanical property is compression energy (y3). 

15 19. A method according to any one of claims 15 to 18 wherein the calibration 
equation is determined from laboratory data establishing a correlation 
between reflectance and the biomechanical property, 

20. A method according to any one of claims 1 to 19 wherein an additional 
property of the feed is also determined. 

20 21 . A method according to claim 20 wherein the additional property of the feed 
is digestibility of dry matter in vivo or in vitro. 

22. A method for determining feed quality, the method comprising the steps of: 

(a) subjecting the feed to infrared radiation to obtain spectral data; 
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(b) using the spectral data to determine a biomechanical property of 
the feed; and 

(c) using the biomechanical property obtained in step (b) to determine 
feed quality; 

5 whereby, the biomechanical property of the feed and thus feed quality is 

determined on the basis of the bond energies of the chemical constituents 
of the feed. 

23. A method according to claim 22 wherein the feed quality is determined as a 
measure of voluntary feed consumption (VFC). 

10 24. A method according to claim 22 wherein the feed quality is determined as a 
measure of forage consumption constraint (FCC). 

25. A -method -substantially as herein described with reference to the 

description of the examples. 

26. A spectrometer configured to carry out the method according to any one of 
15 claims 1 to 21 wherein the spectrometer is adapted to receive a sample of 

feed and determine a biomechanical property of the feed. 

27. A spectrometer configured to carry out the method according to any one of 
claims 22 to 24 wherein the spectrometer is adapted to receive a sample of 
feed and determine the quality of the feed. 

20 28. A spectrometer according to claim 26 or 27 further comprising a data 
processing means for determining the biomechanical property or the quality 
of the feed. 
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Figure 2a 
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Figure 2a (cont'd) 



Energy required to comminute C^J/kg DM) 




2^ 



•jQ^ Digestibility of dry matter in vivo (%) 



S 



40 



Stepwise 2,5.5 



SS 
Actual 



SUBSTITUTE SHEET (RULE 26) 



A1> 



wo 97/21091 



PCT/AU96/00776 



4/11 



Figure 2a (cont'd) 
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Figure 2b 
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Figure 2b (cont'd) 
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Figure 2c 
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Figure 2c (cont'd) 
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Figure 2c (cont'd) 
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